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PREFACE. 

TJ1HE lectures contained in this booklet wen* all delivered 
practically extempore. That they see the light of day, is 
chiefly due to the kindness of Mr. Corakli Prasad, M.Sc , and 
Mr. iBadrinath Prasad, M.Sc., Assistant Professors of Mathematics 
in the Benares Hindu University, of whom the former took down 
notes of the last three lectures and the latter of the first three. 
The copy for the press was prepared by me from these notes 
with the help of Mr. Mewa Ram Saksena., M.A., B.Sc., of the 
Shyamsundar Memorial Intermediate College. Ch and ansi. Dr. 
Bibhutibhusan Datta, D.Sc., Acting Head of the Department of 
Applied Mathematics in the University College* of Science, Cal- 
cutta, kindly helped me in correcting the proof sheets. 1 take 
this opportunity to express my thanks to all these friends and 
former pupils. 


April. 1924. 


CANES FI PRASAD. 
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FIRST LECTURE. 


[Delivered on fche 10th January, 1921.] 

INTRODUCTION. 

The chief reason for my choosing “ The place of partial 
differential equations in Mathematical Physics” as the subject for 
these lectures is my wish to inspire in my audience a love for 
Mathematics. Before entering into details, however, I shall give a 
brief historical account of the application of Mathematics to natural 
phenomena. 

Mathematicians devote themselves to their subject chiefly with 
one of the two objects : either to study natural phenomena in con- 
nection with matter, space or time ; oi for purely aesthetic pleasure. 
A good il lustration of the latter is afforded by their attempt to find 
all the prime numbers less than a certain given number. In ancient 
times people studied Mathematics mainly because of some mystery 



which they attached to it and not for its application in studying 
natural phenomena ; in this connection, we may note, for example, 
the doctrines of Pythagoras or the Sulvasutras of Apastamba. 
Among the earliest writers who devoted themselves to Mathematics, 
were Archytas of Taren turn (428-327 B.C.), Aristotle (3 84-3 2 2 B.C.) 
and Archimedes (287-212 B.C.). Archytas applied his knowledge 
of Mathematics in building up the science of Mechanics. After 
Archytas, Aristotle did something in this respect. But he had 
strange ideas about some elementary truths ; for example, he 
supposed that for a falling body the rapidity of its motion depend- 
ed upon its weight, according to him a heavy body falling more 
rapidly than a lighter one. Another erroneous idea which Aristotle 
had was that if a stone were thrown up at a certain angle, it 
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would go first upwards in a straight line at that angle and then 
fall vertically downwards, the path of the projectile consisting of 
two straight lines, as in the accompanying figure. Archimedes 
discovered the principle of the lever and the conditions of the 
floatation of a solid in a liquid. Galilei (1564-1642) was the first 
to acquire clear ideas about the first two laws of motion. He also 
stated and proved that the path of a projectile was not a pair of 
straight lines but a parabola. After Galilei, Mechanics became 
gradually a science. A large number of results were discovered 
between 1642 and 1700. 

In the beginning of the 18th century a new epoch was created ; 
for, people began to think of natural phenomena other than 
ordinary motions, such as for instance the state of a hot body. 
This gave rise to the use of partial differential equations. Let there 
be a solid body which conducts heat. The amount of heat passing 


Y 



across a sheet at P perpendicular to the r-axis is required to be 
found out first. The amount of heat conducted will be propor- 
tional to the gradient of the temperature and is equal to 


— k 


9V 
dx 5 


which is the first germ of truth in the Analytical Theory of Heat 
and was discovered about the year 1747; the partial differential 
equation governing the state of a hot body was discovered much 
later. Another important discovery was made in connection with 
the study of sound. People had already acquired some idea of 
harmony, and the formula of Taylor, 


viz. 
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was roughly known to Aristotle and the anoient Rishis of India. 
In 1636 Mersenne gave the above formula as the result of his 
experiments. But all such results were at first more or less of an 
empirical nature. The equation arising in the problem of the 
vibrating string, viz. 

d*y T d l y 
dt 2 p dx l 

was given by D’Alembert as late as in 1747 and was the first partial 
differential equation to be carefully studied. It was solved in 
T 

two ways. Putting — = a 2 , 

y = f(x + at) + F(x - at ) , 

was the first solution given by Euler and D’ Alembert. Another, 
which was given by Bernoulli, is 


y=ZA n 


nirx 


sjn 


l 


cos 


n-rrat 
~ ' 


fn appearance, Bernoulli’s solution is periodic whereas tho 
former solution is non-periodic. Therefore the two solutions 
formed the subject of a most interesting controversy. Fourier 
showed that Bernoulli’s solution was also right and that, contrary 
to the views of Euler, it can be used for the case of an arbitrary 
displacement. 

It is the deliberate opinion of modern authorities that before 
the use of partial differential equations, there was no sound basis 
upon which the science of Physics could rest. Soon after the time 
when partial differential equations came to be applied to the study 
of natural phenomena, Fourier made the following statement 
which will remain true for ever : — 

“ The analytical equations, unknown to the ancient geometers, 
which Descartes was the first to introduce into the study of curves 
and surfaces, are not restricted to the properties of figures, and to 
those properties which are the object of rational Mechanics; they 
extend to all general phenomena. There cannot be a language 
more universal and more simple, more free from errors and from 
obscurities, that is to say more worthy to express the invariable 
relations of natural things. 

“ Considered from this point of view, mathematical analysis is 
as extensive as nature itself ; it defines all perceptible relations, 
measures times, spaces, forces, temperatures ; this difficult science 
is formed slowly, but it preserves every principle which it has once 
acquired ; it grows and strengthens itself incessantly in the midst 
of the many variations and errors of the human mind. 

“ Its chief attribute is clearness ; it has no marks to express 
•confused notions. It brines together phenomena the most diverse. 
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and discovers the hidden analogies which unite them. If matter 
esoapes us, as that of air and light, by its extreme tenuity, if bodies 
are placed far from us in the immensity of space, if man wishes to 
know the aspect of the heavens at successive epochs separated by 
a great number of centuries, if the actions of gravity and heat are 
exerted in the interior of the earth at depths which will be always 
inaccessible, mathematical analysis can yet lay hold of the laws of 
the phenomena. It makes them present and measurable, and 
seems to be a faculty of the human mind destined to supplement 
the shortness of life and the imperfection of the senses ; and what 
is still more remarkable, it follows the same course in the study of 
all phenomena ; it interprets them by the same language, as if to 
attest the unity and simplicity of the plan of the Universe, and to 
make still more evident that unchangeable order which presides 
over all natural causes.” 

We have the following important partial differential equations 
that occur in the study of natural phenomena : — 


( 1 ) 


<9*r 

IF 


S 2 V 

<9.H 


( 2 ) 

(3) 

(4) 

(5) 


(9 V 


d 1 1 

d* V 

d* F, 

dt ~ 

c* ( 

dr 1 + 

« y 1 

dz 2 1 

d‘ l V 

- / 

d 2 v 

d‘V 

<9 2 F 

dt* ~ 

' ( 

Ox* + 

W 

+ dz* 

d*V 

c7* V 

d*V 

() 


fix* + 

~dy* 

+ a7 

— V. 



(VV 8*F cl' l V WT7 

+ + + x z V 

ftr* dif dz* 


r) 4 F $iy 

• + — — 47Tp, 

cty 2 dz 

<9 2 F 8F_fl 2 
fl/*“ + H St “ dx* 


As an illustration of the development of the rigorous method 
for the study of natural phenomena, I will sketch an outline of 
the history of the proof of Poisson’s equation. 

^? 4 F = — 47 r/>. 


The proof of Poisson’s equation, usually given in ordinary 
text-books, is not quite correct. The proof is based on enclosing 
the point (x, ?/, z) in a small sphere where the density p is taken to 
be constant. F = F, -j- F* where F 2 is the potential due to the small 
sphere and F, is the potential due to the rest of the attracting 
body. and thus v*F = r a F*. Now if you postulate that p 

is constant everywhere in the small sphere, then of course Poisson 
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proof can hold true. But in case p is variable, this proof will no 
longer hold. Green took up this question in 1828 and said that p 
is not constant but its variation is very small in the sphere and is 
therefore negligible. The first rigorous proof was published in 
1§40 by Gauss who took c and its first partial differential coeffi- 
cients to be continuous. Now suppose F, the potential at P(.r, y, z), 
is due to a solid for which p at any point Q is 

cos a 0 

log- 

r 

where PQ = r , the radius vector and 0 is the angle made by the 
radius vector with the axis of x. Tn this case Poisson’s equation 
is not valid at P ; for, 


d*F 

dx %Z 


d*F 
dy * 


d l V 

dz 2 


has no meaning. I invite your attention to Professor Petrini’s 
generalization of Poisson’s equation, viz. 


a fl 
=o, L ^ 


Lt 
h 

k= 0, 
1 = 0, 


1 K dV(x + h , y , z) dF j 

h I 


1 V BV(x y y + k y z) dF 


dx 


dad 


+ i\ 


1 \ BV{x,y,z + l\ dV 
l l dz dz 


dy 

j- J = - 4ir/>, 


_ I 

dy i 


where h y k y and l tend to zero without the ratio of any two of 
them tending to zero or infinity. There are various advantages of 
this generalization, but in my paper 1 “ On the failure of Poisson’s 
equation and of Petrini’s generalization/’ I indicated a case in 
which this generalization fails. Take for example 

/> = cos log 


then there is no limit in this case except for very special values of 
h, k , l y and the generalization as well as Poisson’s equation fails. 

The application of the modern Theory of Functions of a real 
variable makes the use of partial differential equations more logical. 
Also the application of the Theory of Integral Equations is very 
useful. A typical integral equation is 

b 

<j)(x)=f(x) + X { K(x, 


* Bulletin of the Calcutta Mathematical Society , Vol. 8, 1918. 
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where <f> is the unknown function to be determined. The first im- 
portant paper which laid the foundation of the Theory of Integral 
Equations as a systematic branch of Mathematics was published 
by Professor Ivar Fredholm in 1900 and was later published in the 
Acta Mathematica in an extended form in 1903. Integro-differential 
equations of a somewhat complicated type are useful in the “ Me- 
chanics of Heredity,” i.e. those branches of Mechanics in which the 
state at a particular point of time is dependent not only on the 
immediate past but also on the whole period antecedent to that 
point of time; as for example in the phenomena of Hysteresis. 
Integro-differential equations are equations which involve not only 
the unknown functions under signs of integration but also the 
unknown functions themselves and their derivates. A paper on 
the subject of integro-differential equations by Professor Volterra 
appeared in the Acta Mathematica , Vol. 35. The subject of integral 
equations is of essential importance in studying the phenomena of 



Yj 

radiation. The law of Kirchoff, - = constant, cannot be proved 

a 

without the help of integral equations. The first rigorous proof 
was given by Professor David Hilbert in 1913. The kinetic theory 
of gases cannot be rigorously studied without the help of integral 
equations (see Hilbert’s paper in the Mathematische Annalen , Vol. 
72, 1912). 

There can be no better instrument for studying natural 
phenomena than integral equations, partial and ordinary, and 
integro-differential equations. The truth of this remark can be 
easily seen by considering the problem of Dirichlet : To find a 
solution of v # F = 0 inside a given region (Fig. 3) when, on the 
surface bounding the region, F = / (as, y , z), a prescribed function. 
By the help of the integral equation 

+ p {P) *. 
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the quantity p can be found as the strength of a double layer whose 
potential is the required solution of the problem. Integral equa- 
tions are also most serviceable in studying the problem of the 
conduction of heat in a body. The well-known partial differential 
equation for the linear conduction of heat is 

dF_ % d*V 
St C Sx % ' 

But this equation is invalid for some very simple cases and has to 
make way for a partial integral equation. (See G. Prasad’s mom oh* 
“ Constitution of matter and analytical theories of heat ” in the 
Abhandlungen der Koniglichen Gesellschaft der Wissenschaften zu 
Gottingen , 1903). 
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THE PARTIAL DIFFERENTIAL EQUATIONS OF VIBRATORY 

PHENOMENA. 

The most important partial differential equation met with in 
the treatment of ordinary vibratory phenomena is of the form 

B 2 V d*V 

Bt 2 \ Bx 1 + By 1 Bz’ ) 


Vibratory phenomena of electrodynamic origin are governed 
by the equations of electrodynamics : — 


= ?y 

By B ; 



Sy 

dx 5 


4 7 TW = — 
ox 

B a BZ 
Bt By 

Bb BX 


da 

dy' 
BY 
Bz 9 

BZ 


Bt Bz Bx 


__Bc Y BX 

Bt Bx By 9 


where the components of magnetic force at the point ( x , y , z) are 
denoted by <*, y, the components of magnetic induction by a, b , c, 
the components of electric force by X, F, Z, and the components 
of total current strength by u, v, w . These six equations represent, 
in symbols, the laws of Faraday, and their solution reduces to the 
solution of a single partial differential equation of the type 


B*F 

bF 


= c % 


/B l F t B Z F B*F\ 
\Bx 2 By 2 + Bz 2 ) 


+ 4 -nchi . 


In the study of many kinds of elastic vibrations, we come 
across a partial differential equation of a similar type. In the case 
of the vibrations of an elastic plate, the motion is given by solving 
an equation of the form 
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d^w 

IF 


V\^w) , 


where, as usual, stands for 


il 

fix 1 + fit/ 1 + fiz 1 ‘ 


I proceed to give you two examples to show how partial 
differential equations are utilized to determine the vibrations in 
concrete cases : — 


First consider the motion of air in a rigid spherical envelope. 
Transforming the equation 


f) z V 

~ = c 2 r 2 F 
W 


into polar co-ordinates, we have 

»y 1 S iY,. in I 

d t % r % ( Sr\ fir/ sin 6 d0\ dO / sin 8 # d<p 2 » 


Putting V=U (r, 0 , <f>) . cos kct, we find that the equation is 
satisfied by 

Jm i&r) 

<p), 

v r 


where $„($, </>) is any surface harmonic of degree n. Certain bound- 
ary conditions are to be satisfied in every case. Here 

J n + l 2 (kct) — ^ 

>/ a 

a being the radius of the envelope. This gives the suitable values 
of k and thus we have an infinite number of solutions each satis- 
fying the differential equation and the condition at the boundary. 
The sum of such solutions 


J n + Ak A r) 

V = 2 ( 0 , ?) cos k ct 

V> A 

must also satisfy the differential equation and the boundary con- 
dition and gives the motion, the constants in 8„ tK being so chosen 
as to give the prescribed initial state when t is equated to zero. 

As a second example, consider the vibrations of a circular 
membrane. 

,/ar d 2 v\ 

df l 2 \]dx 2 + Sy 2 / 

is the equation of motion. On transforming it into polar co-or- 
dinates, we have 

d 2 V 


B 2 V 1 dV 1 d*V 
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Proceeding as in the first example, 


SHJ 1 SU 1 
Sr 2 + r Sr + r 2 


+ k*U-=i). 


The typical solution of this equation is 

U = J* (&r) (.4 U cos (9 + sin n 0 ). 


Since at the boundary r = a, there is no displacement, U is 
zero there and therefore k is a root of the equation 

J»(ka)= 0. 

The motion is given by 

V = 2 J n (k K r) cos k A ct (A A cos nO +• A sin nd ) , 

», A ' 

the constants being determined by the prescribed initial state. 

Now I want to give you some criticisms, from the point of 
view of the theory of functions of real variables, on the usual 
method adopted for the study of vibratory phenomena. $or the 
sake of simplicity and fixity of ideas, I will consider the simplest 
vibratory phenomenon, viz. the vibrations of a string. 

The partial differential equation of the motion is 

St* Sx * ’ 

T ... 

where c* stands for — and may be taken to be unity if the units 

are suitably chosen. Let the string be fixed at its extremities^ 
.t=o, #=7r and let it be set in motion by a displacement y=f (x). 
Then the motion at any subsequent time t is given by 

X 

y=. 2 .4„sin nx cos nt , 

n— ] 

the constants being chosen in such a way as to satisfy the initial 
condition . Therefore 

X 

f(x) = 2 ^4„sin nx 


and 


7T 


fix) sin nxdx. 


The flaws in the above method for determining the motion of 
the vibrating string were pointed out by me in a paper 1 “ On the 


t Se© Vol. 6 of the 4< Bulletin of the Calcutta Mathematical Society 1916. 
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vibrating string with an infinite number of edges.” The following 
are the two flaws : — (1) We always assume the existence of 


dt* 


and 


S*y 

dx % 


fiffc every point. (2) It is taken for granted that if each term of an 
infinite series satisfies the differentia] equation then the sum of the 
series will also satisfy the equation. 

Both the assumptions are unjustifiable. For, a string may 
have an edge so that at that point 

dy d*y 

— , what to speak of — , 


will be non-existent ; also what is true of a finite number of terms 
may not be true of an infinite number of terms. These objections 
I proceed to elucidate at some length. 

Consider the string which is under constant tension I 7 . Ac- 
cording to Newton’s second law, study the change of momentum 
of a small element of the string in the period of time 

(t\ - t , t\ + T )- 

Let the ^-co-ordinates of the ends of the element be a^-A, aq + A. 
Then the change in its momentum is 



m i - h 




It must be balanced by the change in the force which is 


t x + r 



+ h 


t x -r 



U) 


Only the assertion of this equality gives the true equation of the 
phenomenon as an integro-differential equation. 

What is usually done is to consider 2 A as getting smaller and 
smaller and also 2 t getting smaller and smaller. Then the appli- 
cation of Rolle’s theorem shows that the integro-differential equa- 
tion is equivalent to 


2A . 2r . 



t^ti + 6' r, 

+• $ih 


= 2t . 2h 


Z(*v\ 

P 


+ *2 t > 
x=*xi+ Qgh 


inhere the 0’s are all proper fraotions. 
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The conclusion, that 



is not always right ; if there is discontinuity then the above equa- 
tion may be invalid. 

Now I proceed to elucidate the second flaw by taking the 
initial state of the string to be y = f(x), where 

f(x)-x iorQ£x£~ , and / {x)=n-x for ~£x£tt . 

The ordinary procedure gives the displacement y at any time as 

4 ^ { _ ] sin (2n — l)x cos (2n- 1)£ 

~ ' (2w-l)* ’ 

the assumption being that 

57 ~= — - 2 ( — 1 )” -1 sin (2n — 1)# cos (2 n— 1)/, 

it n= , . 


4 * 

and x~ = — 2 (-l) tt 1 sin (2n — l).r cos (2n— l)t 

(jX / 7T w — 1 

and therefore the equation 

d 2 y 8Hy 
dt 2 ~~ dx 1 

is satisfied. But the assumptions are unjustified because the series » 
to which 


S l y <9 2 y 

to as well as — , 


is equated has no meaning , it being not convergent . Thus both 


d' 2 y _ 8 2 y 

& r* and w 


have to be investigated by a different procedure and their equality 
established by a different method. 

How to study then the phenomenon of the vibrating string ? 
This question was answered correctly by the late Dr. A. Harnack 1 
in 1887. If 

8 l y , 8 l y 

y <*.*>., yi and w 


1 See his paper, “Ueber die mit Ecken behafteten Schwingungen gesoannter 

Saiten” ( Mathematische Annalen, Bd. 29). 
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are finite and continuous then it was proved by Harnack that the 
most general motion of the string is given by 


Proof 


?/= 2 A„ sin nx cos nt. 

n-1 


Since is existent and finite, according to Fourier’s theorem 


y (x, t) must be equal to 


« w (/) sin nx , 


=-. I y (x, t) sin nx dx. 


where 


[Here it may be noted that every continuous function is not 
expansible in Fourier’s series. This was first shown by Du Bois 
Reymqnd who gave an example of a continuous function the 
Fourier’s series corresponding to which is not convergent.] 
Integrating by parts 


sin nx dx — —n* \ y sin nx dx = — n % - a n (t). 


Since — is continuous, we have 




( ay . 

sin nx dx— 1 ~ sin nx dx , 


which equality is true because ~ is continuous. Therefore similarly 

tit 


V ( , . 

w \ y s * n nx dx = \§^ sl 


sin nx dx ; 


and this latter integral equals 


sin nx dx , 


because 


d 2 v _ d*y 
dt 2 dx 2 


Hence 
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Therefore a n = A n cos nt + B n sin nt. 

If the string is simply displaced initially, then from 

QO 

?/ = 2 (A n cos nt + B n sin nt) sin nx 

n = o 

we find that the B’s must be all zero. 

Therefore the general solution is 

x> 

y= 2 A n cos nt sin nx. 
i 


1 proceed to consider the case when there is an edge at a point 
P on the string so that at that point ~ is non-existent and conse- 
quently is also non-existent there. 

OX 


Denote the ordinates by y ]t y 2 in the portion APB and by 
y' v y\ in the portion AP'B. The usual phenomenon of wave mo- 
tion is that every wave propagates itself. Consider the small ele- 
ment of the string round P during a small period of time r. This 



dc 

element is — * where c is the abscissa of P ; the element undergoes 
dt 

a change 

dc 

^ dt T \ dt dt y 

in momentum and the change is brought about by the impulse 

T J§M&}. 

( OX OX J 

T 

acting for the period of time r. Therefore, since - = 1 , the units 

being suitably chosen, at P 

dc 5 dy, dy t > _dy 2 dy ] 
dt l dt dt » d% dx 


( 2 ) 
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As the string is continuous y l (c, t) = y 2 (c, t ) ; therefore 

y i dc 8y ] dy 2 dc dy% 

3c ' + dt dc it dt 


This gives 


dy, dy t _(k/dy i dy l \ 
dt dt dt\ dx dx ) 


( 3 ) 


In order that the two equations (2) and (3) be consistent, — 

at 

must be equal to 1. Therefore the equation, which must be sa- 
tisfied at the edge in the place of the partial differential equation, 
which has no meaning there, is 



The question is now whether in the present case the general 
solution holds ; it will be easily seen that the required condition 
at the edge is satisfied by the solution and therefore the answer 
to the question is in the affirmative. 

Cases of strings with three or four or any finite number of 
edges have been considered by Harnack. But the case of an in- 
finite number of edges was treated for the first time by me in the 
paper already referred to. As a simple example, it is proved there 
that 

2 « C n 

y = - z — r cos nt sin nx 
it 

gives the motion of the string, vibrating with the infinite number 
nf edges, the abscissae of which initially formed the set 


G n being equal to 


(-D 



2 


ac 

2 



Utx 


2 ' 1 


*7T°r(-l) % 


1 M7T 



according as n is even or odd. 


I will conclude this lecture by showing how the theory of 
integral equations may be applied in the study of the problem of 
the motion of a vibrating string. 

Let us consider the case of a string of variable density. Then 
in the equation 


P (x) 


^y =T y 

at 1 dx * 


(4) 


p is not constant. Lord Rayleigh could not integrate this equation 
exoept for very special types of p. 
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The general solution can be obtained by the help of an integral 
equation. For, put 


2 / = cos nt . U(x). 


Then the equation (4) becomes 


<9 2 £7 n 2 
to* + T r(x) 


U = 0 , 


and it can be shown without difficulty that this last equation to- 
gether with the boundary conditions 

£7(0) = U(tt) = 0 

is equivalent to the integral equation 

7T 

0 

where G ($. x) — — for 

7r 

G (£ x) = — — for 

7T 
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THE PARTIAL DIFFERENTIAL EQUATIONS OF THE ANA- 
LYTICAL THEORY OF THE CONDUCTION OF HEAT. 

After the discovery of the experimental foundations of the 
elementary laws of conduction and radiation by Newton, Lambert 
and Biot, a good deal of time elapsed before the discovery by 
Fourier of the following differentia] equation as representing the 
phenomenon of the conduction of heat in an isotropic solid : — 

dV 2 (d*V d z V_ (PF\ 

(9 f ° \ dx 1 + + Sz % ) ’ 

where *c 2 is the coefficient of diffusivity. A number of special 
problems were solved by Fourier and his contemporary Poisson. 
For example, take the problem of the cooling of a solid sphere, 
whose surface is maintained at a constant temperature, say zero : 
the initial state being given by 

V(x, y, z , 0 ) = f(x, y, z) . 

The particular case in which f(x, y , z) is a function of r alone 
was solved by Fourier ; the general case was considered by 
Poisson. To determine the state inside the sphere at time t , we 
.begin by transforming the partial differential equation into polar 
co-ordinates. Thus 


av a ri ( d ( .av\ i d ( . , av\ 
1F =C b* I SP V aP) + ^o Tv V sin 9 w) 

Putting 

V = e ' c * m U (r, 9, <p) , 


_ &v 

e a<p 1 


i] 


we have 

d / 2 SU\ 1 d / . _ 

I r l — I + - — ^ Tn ( sm 0 
or \ dr / sin 0 <90 V 

of which the typical solution is 


SJJ\ 

&e) 


i dHJ 
+ sin % 6 d<p 8 * 


r*k*U = Q 


+ ± (hr ) 
x/r 




where S n (0, </>) is a spherical surface harmonic of degree n. 
boundary r=a, F is zero ; therefore k must satisfy 

Ju+h (ka) =0 . - 


At the 
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Let the values of k be k v k v k 8 , etc. Then the required solu- 
tion of the problem is given by 


V = 2 

m, w 


Jn 


>h(k»r , 

Vr 


(' M ) 


g A;' 2 w c'H 


the summation extending over all values of n from 0 upwards and 
over all values of m from 1 upwards and the constants in S n , m being 
so chosen as to satisfy the condition relating to the initial state. 

As another example, consider the state of heat of an ellipsoid 
of three unequal axes. Lame in 1837, Heine in 1839 and Mathieu 
in 1870 made important contributions to the advancement of the 
knowledge of the subject. Transforming the partial differential 
equation V 2 F = 0 into ellipsoidal co-ordinates, the equation 
reduces to 

if 

where Re stands for *\/ (a* + 0) ( b ? + 0) (c* + 6) . Taking V td be of 
the form E(^) E(/j) E( v) where E(0) is a function to be determined, 
the equation will be satisfied if 


(fie E(0) = AO + B, 


where A and B are constants. The function E(0) is called Lame’s 
function and the product E( A) E(r) E(v) is called an ellipsoidal har- 
monic. The stationary state of heat of* the ellipsoid can be 
expressed in terms of a series of such products. As regards thd 
non-station ary state of heat, the partial differential equation of 
the phenomenon is 


and this reduces to the form 


where 


V + 


-\) 


(* 0 


+ (*■ 


•«) 




at 


M = 




Mathieu made an unsuccessful attempt to solve the above 
equation. For the special case of an ellipsoid of revolution, Professor 
C. Niven’s results are of some interest. The most successful 
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attempt to solve the problem of the non-stationary state that has 
been made up to now is that of Mr. Bibhutibhuson Datta, who 
obtained, for the general case of the ellipsoid of three unequal 
axes as well as for the particular case of the ellipsoid of revolution, 
approximate results of very great interest and importance, one of 
tfye results being the correction of a mistake of Prof. Niven. (See 
the j Bulletin of the Calcutta Mathematical Society , Vol. 8, and the 
American Journal of Mathematics , Vol. 41). 

T propose now to submit the use of partial differential equa- 
tions to a criticism from the point of view of one who wishes to 
be strictly rigorous. For the sake of simplicity and fixity of ideas, 
I will consider the conduction of heat in an infinite slab bounded 
by the faces x = - tt, x = w, which are impermeable to heat, and I 
will take the initial temperature of the slab to be f(x) where f(x) 
is an even function of x. 

The partial differential equation of the phenomenon is 

SV ffF 
St Sx 1 

where c l stands for the coefficient of diffusivity and may be 
taken to be unity if the units are suitably chosen. Also because 
of the impermeability of the bounding faces, %=. tt, x = n, we 
SV 

have -r- = 0 for x = t as well as for x = ~ Thus the state of the 
ox 

slab is given by 

00 nU> 

V = \a t) + 2 a n e ' cos nx, 0, — ? -n . 

1 

‘where 


9 ( 

a n — ~ \ f(%) cos nx dx, 
0 


The above-mentioned treatment of the problem is the ortho- 
dox one and will be found in every text-book. In Fourier’s 
original papers the statement of the conditions is the same as 
above, with this difference that it is clearly mentioned that 


SV S l V 
St Sx 2 


for every positive value of t including zero and for every value of 
x in the interval ( - tt, tt). 

In order to establish the defective character of the orthodox 
treatment, I take the simple case in which f(x) is given by the 
accompanying graph. 
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At 0 there is no differential coefficient and still by the ortho- 
dox method 

T7 it 4 * cos(2ft— 1)& - {2n~\)H 

If we differentiate the abqve 
series term by term the series 
so obtained is uniformly con- 
vergent as long as t > 0 and 
hence we can differentiate in 
that manner for such values 
of t. Thus for t > 0, 

<9 V (0, t) * ~{2n l)H 

at ~7 e 

Fig. 5. and, making t tend to 0, 

lim d V (0, t) 
t = 0 Si ~°°* 

Now according to the theory of functions of real variables, if 

lim dl'(0 ,/) 
t = 0 St 

is existent, it must be equal to 

lim F(0, /)— F(0,0 ) 

/ = 0 t ‘ 

s v 

Therefore, at the time £ = () the value of at 0 is infinite 

cH t 

and also at that point, there being no differential coefficient of F 
<9*F 

with respect to x , g^-is non-existent ; consequently for x = 0, £= 0, 
the equation 

SV S 2 V 
dt ~~ Sx 2 

has no meaning and the orthodox treatment fails. 

How to study then the phenomenon of the conduction of 
heat ? 

This question was answered correctly for the first time in my 
memoir, “ Constitution of matter and analytical theories of heat ”, 
published by the Royal Society of Sciences of Gottingen in 1903. 
(See Abhandlungen der Koniglichen OeselUchaft der W issenechaften 
zu Gottingen , Bd. 2, No. 1). In the theory, as remodelled by me, 
the equation of the phenomenon is shown to be not a partial 
differential equation but a partial integral equation, and attention 
is drawn to the simple principle that it iR physical continuity 
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which is necessary and not the continuity of the form of an ana- 
lytical expression. I proceed to elucidate these features of the 
theory. 

First , as regards the equation of the phenomenon, this 
follows from the application of the principle of the conservation 
of energy in accordance with which the amount of heat passing 
into a volume must be equal to the amount of heat absorbed to 
raise the temperature of the mass of the volume. In the form 
of symbols, this condition may be written as : — 


/ m 

(A) j dt ' = j [T(X', t)-T(x',0) j dx', 


T(x, t) denoting the temperature. It can be easily shown that 
this integral equation is equivalent to the differential equation 


dT^&T 
dt dx* 

dT d*T 

if — and are finite and continuous. For, by the mean value 

theorem, 
h 


ft 




BXy 


dxt 


0, and 0 2 being proper fractions. Similarly 


+ 0J (a? a — ap] ) 
t ~ ^1 t\) 





From (A), the two integrals given above are equal; consequently, 
proceeding to the limit by making t 2 - ^ and x. z - x l tend to zero, 
we have 

d_T_d*T 
dt ~ dx % 


for x=zx lf t=t l9 if these differential coefficients are finite and 
continuous. 

Secondly, T(x , 0) need not be derivable from the analytical 
expression, which equals T(x , t) when t>0, by putting f=0 in that 
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expression. Thus, in the remodelled theory the convergence of 
the Fourier’s series 

QO 

ha () + 2 o n cos nx 
1 

has no place. On the other hand the continuity of T(x, t ), i.fe., 
the equality of 

lim V(x, t) to f(x), 

1 = +0 


follows at once from the continuity of f(x) ; and the only restric- 
tion imposed on f{x) i in order that F(:r, t) should be the solution, 
is the following : — 

(B) For each value of x, within the interval (-? r, tt), the upper 


limit of 


SV 

Sx 


in the neighbourhood of £ = 0is finite. 

SV . 


O V 

The necessity for the finiteness of — is obvious ; for, since 


dV 

-k — is the amount of heat which flows across an unit area of the 
dx 


plane at x in unit time, if this amount becomes infinite it will 
be absurd. Now for t > 0, the series for V(x 9 t) being uniformly 
convergent and the series obtained by differentiating it being also 


SV . 


uniformly convergent, ^ is finite. 


Hence we have to consider only 


the neighbourhood of £ = 0. 

The equation (A) is satisfied by 

T(x , t) = V(x, t), t > 0 

as long as f(x) is finite and integrable. For, for i> 0, the series 


SV 

for and V can be integrated term by term, the series concern- 


ed being all uniformly convergent. 

SV 

As regards (B), it can be satisfied even if -gj- makes an infinite 
number of oscillations as t tends to 0 ; for example, if 


S_V 

St 


X 



makes an infinite number of oscillations with indefinitely 


increasing amplitudes as t tends to zero. Thus theoretically the 
slab may, at certain points, become hotter and cooler an infinite 
number of times. 
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I conclude this lecture by pointing out that a number of in- 
teresting problems in the theory of the conduction of heat can be 
best solved by the use of the theory of linear integral equations. 
For example, one of the earliest applications of that theory is by 
Professor Picard, to the problem of the cooling of a ring. The 
partial differential equation of the phenomenon is 


<9F 

dl 


d 1 V 

c 2 -FT- h * V ‘ 
dx* 


This differentia] equation can be solved in terms of the solu- 
tions of an ordinary differential equation of the second order which, 
in turn, is equivalent to an integral equation of the type 

U(i) = F(t) + U(x) A Lx) G(x, Z)dx , 
where TJ is the unknown function. 



FOURTH LECTURE. 

[Delivered on the 3rd February, 1921.] 


THE PARTIAL DIFFERENT TAL EQUATIONS OF THE THEORY 
OF GRAVITATIONAL ATTRACTIONS AND THOSE OF 
ELECTROSTATICS AND MAGNETOSTATICS. 

In the three lectures which I delivered last month, I first gave 
a historical review of the present theory of the partial differential 
equations of Mathematical Physics together with some criticisms 
and then dealt in succession with the equations of vibratory 
phenomena and of the analytical theory of the conduction of heat. 

I wish to preface the three lectures which T intend to deliver 
in this month by explaining to you the two kinds of Mathematical 
Physics. This cannot be better done than in the following lyords 
addressed by Professor Vito Volterra. of Rome at the inauguration 
of the Rice Institute of America in 1915 : — There are two kinds 
of Mathematical Physics. Through ancient habit we regard them 
as belonging to a single branch and generally teach them in the 
same courses, but their natures are quite different. In most cases 
the people who are greatly interested in one despise somewhat the 
other. The first kind consists in a difficult and subtle analysis 
connected with physical questions. Its scope is to solve in a 
complete and exact manner the problems which it presents to us. 
It endeavours also to demonstrate by rigorous methods statements 
which are fundamental, from mathematical and logical points of 
view. 

“ I believe 1 do not err when I say that many physicists look 
upon this mathematical flora as a collection of parasitic plants 
grown to the great tree of natural philosophy. But this disdain 
is not justified. In the evolution of Mathematical Physics these 
researches are to play an ever increasing part 

“ The other kind of Mathematical Physics has a less analytical 
character, but forms a subject insuperable from any consideration 
of phenomena. We could expect no progress in their study 
without the aid which this brings to them. Could any one imagine 
the electro-magnetic theory of light, the experiments of Hertz and 
wireless telegraphy, without the mathematical analysis of Maxwell 
which was responsible for their birth ? ” 

To-day I shall confine myself to the first kind of Mathematical 
Physics and the lecture will be mainly destructive in character. 
But to-morrow I shall speak on the constructive side of things, 
when I shall take up the second kind of Mathematical Physics. 
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The following are the most important equations of gravita- 
tional attractions, Electrostatics and Magnetostatics with the 
respective years of their discovery : — 

(Laplace’s equation, 1781) 

(Equation of Coulomb and 
Poisson. 1811) 


(Poisson’s equation, 1813) 


(Helmholtz’s equation, 1853) 

a being the strength. 

[Xhe exact form of the equation (2) was given by Poisson in 
1811. The last equation gives the potential difference at the 
surface of a double layer of strength «-.] 

Let us take up Poisson’s equation 

WV d z V »V 

&*S?*aF = -* wr - 

The proof as given in text-books postulates /> to be a constant 
throughout a sufficiently small sphere round (x, y, z). Let u and 
U denote respectively the potential due to the small sphere 
round (x, y, z) and that due to the remaining mass. Then 

V=U + u. 

By Laplace’s equation X' l U — 0 

and consequently r 2 P= r 2 u. 

Taking p to be constant 

u =5 2np (a 2 — J r 4 ) 

at a point (£, v, £) in the sphere distant r from ( x ; ?/, z), a being the 
radius of the sphere. Differentiating two times with respect to 
£, y 9 t we have at the centre 

V l u = — (2 + 2 + 2) = — 47 rp. 

3 

In 1828, Green gave a proof more or less like the above. The 
first really rigorous proof was given by Gauss in 1840. Gauss’s 
proof, which I will reproduce presently, is based on the supposition 
that p is finite and continuous together with its three partial 


trv o m v err 
(1) dx 2 + dy* + da* “ ' 
dV dV 


( 2 ) 


f r being the surface density. 

<9^ c^F tPF 
dx 1 + dy‘ l + ctz 1 

p being the volume density. 


= - 4 itp. 


(4) »F +11 -»F-o=-4» V . 
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differential coefficients of the first order. In 1882, Professor Holder 
succeeded in establishing Poisson’s equation on the supposition 
that 

p(t,y,0-p{x,y,z) I <A 


where p and A are positive constants and r is the distance between 
(£, rj }, £) and (a;, y , z). The next investigator to impose still less 
restrictions on p was G. Morera who proved in 1887 that, if 


r 

ft 


p(i- v, 0-p(x- y ■ 2 


} 


dr 

r 


is finite, then Poisson’s equation is valid. In 1891, Kronecker 
attempted to establish the equation with lesser restrictions than 
the preceding investigators. From 1891 to 1900 practically nothing 
was done to improve the proof of Poisson’s equation. In 1900 and 
the next ten years, Dr. Henrik Petrini tried to make the proof as 
little restrictive as possible. Some of my papers deal with the 
criticism of Petrini’s work. Petrini’s chief result is that if 

2-n 1 a 


lira i d<p i (3a- 1 - 1 ) d,- \ — 

A = °J J J r 


exists and is finite, then ^ exists and is finite, p standing for 

cos 6 where <9 is the angle between the radius vector and the axis 

of x : with similar conditions for — and ttt *» fbe additional con- 

dy' 1 oz L 

dition for the validity of Poisson’s equation being the continuity 
of p which condition, however, has not been sufficiently emphasized 
by Petrini. But cases may occur in which p is discontinuous. 
For example, I pointed out in a paper published in the Philo- 
sophical Magazine in 1917 that, in the case in which p is 


cos 


^ log , Poisson’s equation is invalid. Instead of cos ^ log - ) 

we may take, without making any material difference, cos ^ log J ^ 
plus any constant or differentiable function. If, then, you take a 
solid sphere of density 2+ cos ( log -j , say, and suppose it divided 
into thin concentric shells, then, by Newton’s theorem, the attrac- 
tion due to a shell is zero at a point inside it and therefore 
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== 0 at the centre. The mass of the shells between the centre and 
Or 

the point (£, rj, £) is 



4t r ^^2 -f COS ^log t % dt. 

0 

Since the attraction at (£, >7, £) is the same as if the whole of 
this mass were concentrated at the centre, the attraction is 


^ {a +oo 8 (log ^)} *dt 
0 

and is along the radius vector. Now at the centre 


B*u 

TT = lim 
Bf 1 .. t\ 



2 + cos log 



-- lim 
r = 0 


4?r 

T~\/io 



i. + tan 

r 



As r is indefinitely diminished, ~ tends to <» and cos ^ log ~ ^ 

S^u 

makes an infinite number of oscillations. Therefore is non- 
existent at the centre. Moreover the radius vector may be taken 
in any direction. Hence, at the centre, 

8 2 u B 2 u 8 z u 
Ox % ’ By 2 ’ <9z 2 


are all non-existent and consequently the equation of Poisson is 
meaningless. 
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For the sake of its historical importance I proceed to give 
you Gauss’s proof of Poisson’s equation : — 

Let p~f(( 9 v> £)• Then the ^-component of the attraction at 
the point ( x 9 y , z) is 


v C (£-#)/(£, y 9 Z 
X(x, y 9 z) = 1 

J {(t-x)* + (ri-y) % ■ f 


(t- 


z ) 2 } 1 




where dt is the element of volume of the sphere and the integral is 
taken over the whole sphere. 

Also 


X(x + e , y , z) = 




(j-x-e) /(£, y, £) 

{ (f- *- c)» + (!,-»)* + ({- 8)* } J 


( 2 ) 


the integral being taken throughout T, the original volume occupied 
by the sphere. Suppose the sphere to be displaced a little to the 
left so that its centre is at P v i, e. ( - e, 0, 0). Let T { denote the 
volume occupied by the sphere in its new position. Then, as is 
clear from the accompanying figure (7), 



Now, put 
Then 


Fio. 7. 

T-d = T ] -e ] = T iy 


X(x + e, y, «)= 1 

.) {(?-*}» + (.,-*)• + <«- 
T. 


dt 


( 3 ) 
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Also 


X(x + e, y, z) - X(x, y, z) ^ ^ + 


I — i. A- A, 


where Z is the integral in (1) taken over T 0 , A the same integral 
faken over 9, Z, the integral in (3) taken over T n and A, the same 
integral taken over 6 V 

The value of an element of volume in 6 and that of an 
element in 0 j are respectively e cos a and - e cos a dS where a is 
the angle made by the external normal with the x ~ axis. 
Therefore in' the limit, 


y*-*) ffcv, 0- 

*/V 


£-.t 

1* 




taken throughout the surface $ of 7 1 „ ; 


also in the limit 


taken throughout T ( , . 



dfj 

s$ dt 


Therefore 


^Z± d Jld i (£-*)* 

dx dx 2 j r 8 d£ ^ 1 r 4 


and similarly for 






a? * n<i 5? ' 


f)dSy 


Therefore 


Since 



T 0 S 

dt = r 2 sin 9d0d<pdr , 


the volume integral in the above equality will tend to zero as the 
radius of the sphere is made to tend to zero ; and similarly 


f 



will approach -4*/ ( x , y , z). 


Thus the equation 

-4irp 


holds. 

I proceed now to consider the generalization of Poisson’s 
equation given by Dr. Petrini. Dr. Petrini was the first to note 
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clearly that the continuity of p was not sufficient for the validity of 

cos ^ 0 

Poisson’s equation. Taking — , he showed 1 that Poisson’s 

d*F d l V d 1 V 

equation fails because is + oo while ^ and are both equfil 
to - oo . Dr. Petrini was thus led to the following generalization : — 


lim r 1 

*i=°. I T 

h Q =o, L 


£dF ( x + h u y , z) d V (x, y , z) 


1 fdV (x, y + z) 


dx 


dx 


K 


l 


dy 


h .=0 


(IV ( x , y, z) \ 1_ 

J + h,_\ 


1 f dV (a:, y. z + h~) dV(x,y,z) 


}]- 


r hj J ' h, h d z dz 11 1 ’ 

provided that while h { , h % , h, are made to tend to zero in any 

manner. ^ , 7- 1 , r should tend neither to zero nor to infinity. Dr. 

' h 9 h % hv 

* cos ^ () 

Petrini proved the validity of his generalization for /<=-— — p and, 

log - 
r 


indeed, for every case where /> is continuous. But he did not con- 
sider the cases in which p has a discontinuity of the second kind. 
Such cases have been studied by me in the paper referred to above. 
I have proved, for example, that when 


1 du 
r dr 




which shows that the limit in the generalization of Petrini will 
exist and be equal to the limit of 

J cos ( log + tan + cos (log ~ + tan ^ ) 

^/io l v "1 37 v K 37 


+ cos 


f 1 -I 

log ~ + tan 

1 h B 



if the latter limit exists. Now put a A,=A 2 and / 3 A i *=A 3 . Then 
1 - 1 1 

denoting log + tan - by log H , the expression within the crook- 
ed brackets becomes 


1 For another method than Petrini’g, see my paper in the Phil . Mag., Vol, 
34 , 1917 . 
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cos (log H) + cos ^log H + log + COS ^log ti + log , 

i. e., cos (log H) | I + cos ^log + cos K)i 

- sin (log H) | sin (log ^ + sin (log | . 

Now the only way in which this expression can have a limit 
is by the vanishing of both 

•^1 + cos (log + cos (log ^ j- and \ sin (log + sin (log j 

which is possible only when a and fj have respectively the forms 
c 2m„ + 2./3 and e 2f».T2*/3 


m and n being integers. Thus, it is proved that the generaliza- 
tion of Petrini does not hold for every discontinuous distribution. 

Another generalization of Poisson’s equation was given by 
Professor C. W. Oseen of Upsala. As I have already taken a 
good deal of your time this evening, I will take up the consider- 
ation of that generalization to-morrow. I will conclude the 
lecture by discussing briefly the equation of Coulomb and Poisson 
and Helmholtz’s equation. 

T have studied the equation of Coulomb and Poisson in three 
papers one of which appeared in the RpMdiconti del Circolo Mate - 
matico di Palermo in 1917, another in the Philosophical Magazine 
in 1918 and the third in the Bulletin of the Calcutta Mathematical 
Society in 1919. If round P, any point on the surface, a small disc 
be described lying on the surface, then for investigating the validity 
of the equation ^ y ^ y 

4 . = — 4t7TCT , 

<9n, 0n 9 

we need only take into consideration the part of V due to the disc. 

Among other results, I have proved that, denoting by r the 
distance of any point of the disc from P, if 

c — cos x(r) when \(r) 

and, further, if 

X (r) ~ log - 
r 

SV dV 

then ^ , — are both non-existent and the equation of Coulomb 
and Poisson fails. 

Results similar to those mentioned above have been obtained 
by me for the validity of Helmholtz’s equation in a paper which 
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appeared in the Proceedings of the Benares Mathematical Society 
for 1920. I have shown, for example, that, if the strength 

(T = cos log ~ , 

then W + 0 and TF_ 0 are hoth non-existent and therefore the equation 
fails. 



FIFTH LECTURE. 


[Delivered on the 4th February, 1921.] 


THE PARTIAL DIFFERENTIAL EQUATIONS OF THE 
THEORY OF GRAVITATIONAL ATTRACTIONS, AND 
THOSE OF ELECTROSTATICS AND MAGNETO- 
STATICS ( Continued ) . 


Yesterday I promised to lay before you the generalization of 
Poisson’s equation given by Professor C. W. Oseen. The general- 
ization 1 is fairly simple and will prove as useful as Poisson’s 
equation. It is this : — 



where dv is an element of the volume a s , ds an element of the 

dV . 

surface S which bounds the volume 12# and ^ is the differential 

coefficient with respect to the normal taken outwards. The only 
condition imposed here on p is that it should be finite and 
integrable. 

It will be remembered that in proving the theorem 
r *V = — 4rrp 


we have to impose more restrictions on p than that it should be 
finite and integrable ; for example, in Gauss’s proof, it is supposed 
that p is continuous and also that p x , p (/ , p r exist and are conti- 
nuous. We see at once how much more simple Professor Oseen ’s 
generalization is. 

Starting from the equation, known in Integral Calculus as 
Green’s theorem, 

one can deduce the modified theorem of Green, viz., 


F(R) = 




1 £F' 
r On 



} 


ds , 


1 See Oseen’s paper, ** Ueber einen Satz von Green und ueber die Definitionen 

von rot und div ’ ’ / Ren.dinnn.ii dal Dirr.nln <m.ntemnt.inn di Pnlermn t. JMl. 101511 
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where // stands for i J (P') and V' for V(P'), P' being the variable 
point (£, y, c), P the fixed point ( x , y , 2), and r the distance PP\ 
This result follows from Green’s theorem which holds when 
v are existent and integrable. To show this, I take U to 

be -pjp . Now yjp becomes infinite at the point P. In order 

to make Green’s theorem applicable, we enclose the point P within 
a small sphere 12 ' of radius a . Inside the region bounded by this 
sphere and S, U and its differential coefficients are finite and 

continuous. Hence, since - satisfies Laplace’s equation. Green’s 

theorem applied to this region gives 



where ds' is an element of the surface S' of the sphere 12'. Now 
let a tend to zero, then, dv being of the order r*dr, the second 
term in the left-hand side of the above equation is negligible. 
Similarly, as ds' is of the same order as r l , the third term in the 
right-hand side is negligible. Also v‘ l ] r '-- 47 r / according to 
Poisson’s equation. Hence we have in the place of (1), 

— 4n j' = 

\l s 

which becomes, in the limit when a tends to zero, 



F(P) = 




i dr 

r dn 


V 


Xh 

dn J 


ds . 


8 

This result is the modified theorem of Green and it says that 
the value of V at any point is the potential due to (1) a distribu- 
tion of volume density p' t (2) a certain distribution on the bound- 
ing surface and (3) a certain double layer on that surface. Thus 
we see that if there is a solid with a known volume distribution 


and prescribed values of V and ^ on the bounding surface, then 

by applying the modified theorem of Green we can at once put 
down its potential V at any point. 

The question arises : Can we deduce the modified theorem of 
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Green by using Oseen’s generalization instead of Poisson’s equa- 
tion ? The question has been answered in the affirmative by Oseen 
who has deduced the modified theorem of Green by applying 
vector analysis and by employing the following theorem: — 

If <p and ^ are two functions, finite and continuous together 
wjth their first differential coefficients, then 


r i* < i( p r 

I (grad <p x grad t) dr - W ^ ds — 471-1 ^ 


The procedure adopted by Oseen is rather complicated, since 
the second differential coefficients of <P may be non-existent, but 


the deduction is an epoch-making result. 

The equation of Laplace, which I have frequently mentioned 
in some of my previous lectures is intimately associated with 
problems in Electrostatics and Magnetostatics. Thus we come 
across this equation when considering the case of a solenoidal 
distribution of electricity and the most important problem associa- 
ted with the equation is Dirichlet’s problem, viz. to find the 
potential at points inside a region when the value of the potential 
is given on the boundary. 

With reference to Diriehlet’s problem, two theorems have 
been formulated and proved under certain general conditions. 
These theorems are : — 


(1) (The Theorem of Uniqueness) : The solution of Diriehlet’s 
problem is unique. 

(2) (The Existence Theorem) : The prescribed value on the 
boundary is arbitrary and corresponding to each value there exists 
a solution. 

That both the theorems hold, was at first guessed from certain 
physical consideration *. For example, in the case of a conductor 
connected with the earth and influenced by external electric 
masses, the value of the potential of the induced distribution on 
the surface is prescribed, being equal in magnitude, and opposite 
in sign, to that due to the masses, and the uniqueness and 
existence of the potential is obvious to a physicist whatever the 
masses may be. But attempts were made to base the two 
theorems on something more rigorous. 

Diriehlet formulated the principle, known as Dirichlet’s 
principle, that, of all the functions, F, which inside the region 
bounded by a closed surface remain continuous together with 
their first differential coefficients and take a prescribed value on 
the surface, there must exist one (or more) for which the integral 

${(£/•(£)*♦(£)> 
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taken over the whole region has a minimum. Accepting this 
principle, we get from the calculus of variations the result that V 
must satisfy Laplace’s equation. 

But Weierstrass pointed out in 1860 that Dirichlet’s principle 
was not always valid and that all that could be said was that tjie 
integral has a lower limit and that it may not have a minimuvn. 
In 1899, Hilbert established the validity of Dirichlet’s principle 
on the supposition that the prescribed boundary value is con- 
tinuous together with the differential coefficients of the first two 
orders. 

Assuming that in a given case the solution of Dirichlet’s 
problem exists and is unique, the question arises of actually 
obtaining this solution. I propose to devote the remaining time 
to-day to the discussion of a number of methods. These are (1) 
the use of Green’s function, (2) the use of fundamental functions 
based on a specially suited co-ordinate system, (3) the method e du 
halayage , (4) the method of arithmetic means of Poincar6 and (5) 
the use of integral equations. t 

Adopting the procedure used in deducing the modified theorem 

of Green from Green’s theorem but taking V to be G + - instead 
of 1 , we obtain for V{P) the simple expression 



here G is regular everywhere inside the region, where it satisfies 
Laplace’s equation, and equals - * at every point on the surface ; 

G + ^ is called Green’s function for the particular region. The 

problem of obtaining the solution of Dirichlet’s problem has been 
thus reduced to that of obtaining Green’s function which may 
be interpreted as the potential due to a unit of electricity at P 
plus the potential due to the distribution induced on S by this 
unit when S is connected with the earth. In some cases Green’s 
function is easily obtained by applying Lord Kelvin’s theory of 
images. For example, let S be a spherical surface of radius 
a. Then the value of Green’s function at any point Q is 

a 1 1 

~ R ' P'Q + PQ ’ 

where R denotes the distance of P from the centre of the sphere. 
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P ' the inverse of P with respect to the sphere ; in fact 5 units at 

H 

P' form the image of a unit at P. 

I proceed now to consider at some length the second method, 
viz. that in which use is made of fundamental functions based on 
a specially selected system of co-ordinates suitable for the parti- 
cular boundary surface. 

As an example of this method, consider the case in which the 
bounding surface is a sphere. Then, transforming v a F = 0 into 
the polar co-ordinates r, 0 , <p, r~a being the bounding spherical 
surface, we have 


<9 / 0 BV\ 1 (9 / . „ HV\ l <9*F „ 

Wr V 


We start by obtaining a typical solution of this equation in 
the form of a product of three factors, one a function of r alone, 
another a function of 0 alone and the third a function of alone. 
Such a solution which is the desired fundamental function, is easily 
found to be 


where 


cos mp 

r n P n m (cos 0) or 

sin mp 


1 fin + in 

(i— 

dfx 


The solution of Dirichlet’s problem is given by 


00 M 

F= 2 1 r n P n m cos + si nmf), 

»*0 »sfi 

where the ,4’s and ZT s are constants which have to be determined 
by the condition that on the boundary, r=a, V takes the prescrib- 
ed value, say /(0, p). Now in actual practice the determination of 
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these constants is not an easj^ affair. Gauss had given 1 a method 
for expanding any arbitrary polynomial (in the co-ordinates of a 
point on the spherical surface) in a series of spherical harmonics. 
Tf <P(x, y , z) be the polynomial of degree n in .r, ?/, z where 


then Gauss put 


x 2 + y z + z* = a*, 




-tin + ti n ^ t + t 4 S„ i + 


ti m standing for a solid spherical harmonic of degree ra, and, by 
operating on the above equation with r 2 , etc., obtained a 
number of simultaneous equations by solving which the S ' s were 
obtained. But no compact and general formula for writing down 
the expansion had been obtained until about twenty years ago I 
gave it in a paper published in the Messenger of Mathematics , Vol. 
30, in 1900. 

The result is r- i f 


{'In — 4/' + 1 ) 


<G.r, y, z) = (- 1) 5 

e o 

where 








2 p ( 2 /*- 2 )... 2 .( 2 w- 2 /x+ 1 ) (I'm- 2 /a- 1 ) . .. 3.1 


A a- a _ 0 — OL — 

dx dy ' ^ bz a (hi 1 + <9/? 2 r)y 2 


The above-mentioned result was obtained by me at first by a 
method different from that of Gauss whose method was applied by 
me to get the result in a paper published in the Mathematische 
Annalen , Vol. 75, in 1912. 

If / (0, <p) is expressed as a function of :r, ?/, 2 , the co-ordinates 
of a point on the boundary, that function can, in its turn! be 
expressed as a sum of polynomials each of which may be expanded 
by my formula, and thus an expansion for /(#, <t>) can be obtained. 
With the help of this expansion, the solution of Dirichlet’s problem 
can be put down at once as was done by me in my paper of 1900. 

As a second example of the method of fundamental functions, 
consider the case in which the boundary is a cylinder. Then, 
transforming V f F = 0 into cylindrical co-ordinates />, s, </>, p — a 
being the boundary, we get 

&V 1 3F 1 &Y cFF 
8? + p dp + ? dp + dz l ~ ' 


The typical solution of the above is 

+ lez cos m p 

e J m ( k P ) or 
sin m0 


1 See his Gesammelte Werke”, Bd. 5, p. 360. 
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this is the desired fundamental function, and the solution of 
Dirichlet’s problem is given by 




„ cos ni(f> + B„ tUt sin mt> 


K kt%z 

c tl , m e + />„, , 


m t 


where the constants, the A' s, ITs, O' s and D’s are determined by 


the boundary condition. 

As the third example of the method of fundamental functions 
consider the case in which the boundary is an ellipsoid. Then, 
transforming r 2 I 7 = Ointo the ellipsoidal co-ordinates A, /*, r, \ = 0 
being the given ellipsoidal boundary 



where 

He = s/ ( a 2 + 0) (//^ -f tf) (6 l + 6/) • 


We suppose, as usual, that a typical solution of the above 
equation can be obtained as a product of three factors, each being 
a function of one co-ordinate alone. Thus we get as the funda- 
mental function 

E( A) E(ft) E(r). 

where 

(lU ~- ; J E(ff) = (mO + n) E(b) . 


which is Lame's equation and can be solved in series The 
required solution of Dirichlet’s problem is given by 

V^A m , n E m%n (\) E m , n ( r ) E m , u (v ) , 

the A " s being constants which are determined by the condition 
that, for A = 0 , V has the prescribed value, say F(/s *). Many 
attempts have been made to obtain a compact and general formula 
for expanding F(g, r) in ellipsoidal harmonics similar to that given 
by me for expansion in spherical harmonics, but so far the attempts 
have proved to be unsuccessful. 

The question of the determination of the most general system 
of three orthogonal surfaces, namely, that of the confocal cyclides, 
for which the fundamental function is the product of three factors, 
as in the cases considered above, was investigated by Professor 
Klein and his pupil, the late Professor Maxime Bocher. (See 
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Bocher’s book, “ Reihenentwicklungen der Potentialtheorie ”, 
Leipzig, 1894.) For the case of an arbitrary bounding surface, 
Poincar6 established in 1896 the existence of fundamental func- 
tions, which, however, are not necessarily expressible as the product 
of three factors, each factor being a function of only one co- 
ordinate. (See his paper in the Acta Mathematica . Vol. 20.) * 

Unlike the previous method, the methode du balayage of 
Poincare 1 uses approximation functions each of which has the 
prescribed value on the bounding surface and each satisfies 
Laplace’s equation more approximately than the one preceding it. 
Briefly, the method may be outlined as follows : — 

Let K be a sphere which encloses the surface 8 and the region 
U s bounded by it and let the prescribed value f(x , y , z) be such 
that / together with its differential coefficients of the first two 
orders is continuous in K : also assume that r a /Z.O in A r , and let <p 
denote 


A 

Then obviously Y *(0 — /) = O in U s . Now construct an enu- 
merable aggregate of spheres K (v K r . . . such that each sphere 
lies in Qs and fills it in such a manner that each point of Q s is 
inside at least one sphere and the spheres become indefinitely 
large in number only in the neighbourhood of 8. Then, if the 
sequence of functions U (v U. 2 , ... is determined by the condi- 
tions that (l) v 2 U 0 = 0 in K n and U 0 equals <p in K - K t) , ( 2 ) 
v a ?/, = 0 in K ] and U { equals U 0 in K - A r ,, (3) r 2 U 2 = 0 in K L and 
V % equals C7, in K -K„ etc., U 0 A U { ^ U 2 . . . and the sequence 
t/ () , U j, U it . . . converges uniformly to a regular potential-function 
in U s : the required solution of Dirichlet’s problem is 

V - ,im + 

The case in which u 2 / >0 can be dealt with by making it 
dependent on the case considered above by suitable superposition. 

The method of arithmetic means as given by Neumann 2 
postulates the convexity of 8 and determines the solution of 
Dirichlet’s problem as the potential of a double layer. If 



1 See hiB paper. “Sur les equations aux d6riv6es partielles de la Physique 
Math6matique ” {American Journal of Mathematics , Vol. 12, 1890). 

2 See his booklet, “Untersuchungen iiber das logarifchmische and New- 
tonsche Potential, ’* Leipzig, 1877. 


itl.r' 
4ir ^ r 


! . < it . 
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s 


then it has been proved by Neumann that the series 
-C + (W-lf 1 ) + (»f r a -lF l 0 + .. . 

s 

converges in ii s , satisfies Laplace’s equation and, with the approach 
of the point P, at which the potential is taken, to the boundary, 
tends to /. Here C is a constant and equals the value of 
Jim W u on the surface. 

n -QO 

The method in which the theory of integral equations is used 
is allied to the above method inasmuch as the solution of Dirichlet’s 
proble/n is obtained as the potential due to a double layer. Let 
the strength of the layer be />, then it is known that 

W e - W,= -4tt p (P) 9 and 

1 cos \ 1/ 

W e + W, = 2 \ P (M) ds , 

8 

where W„ or W \ is the limit to which the potential W of the 
double layer tends according as the point P on the surface is 
approached from the outside or from the inside, P, M, having the 
same meanings as in the figure on p. 6. Adding the two equations, 
we have 

2 W'=- 4 ir P (P) + 2 J p(M) ds . 

s 

But W„ is equal to the prescribed value /(P). Therefore is 
given by the integral equation 

£ 



SIXTH LECTURE. 


[Delivered on the 5th February, 1921.] 

THE PARTIAL DIFFERENTIAL EQUATIONS OF HYDRO- 
DYNAMICS AND THOSE OF ELECTRODYNAMICS 
AND THE THEORY OF ELECTRONS. 


1 propose to speak to-day first on the partial differential 
equations of Hydrodynamics. Less than three hundred years 
ago, Hydrodynamics was not a science at all. In fact the first 
theorem ever discovered in Hydrodynamics is what is known as 
Torricelli’s theorem. This was first stated in the year 1643. 
Torricelli was one of the greatest engineers of his day. He / ound 
from observation that the velocity of water issuing from a small 
hole in a vessel of water is the same as that of a particle falling 
from the level of the water in the vessel, i.e., g 2 = 2 gz, where q is the 
velocity of water flowing out and z is the depth of the hole below 
the water level. The first attempt at a mathematical proof of 
this, was made by Newton. He pointed out that two important 
points were to be kept in mind. ( l) The fluid particles in any 
horizontal section of the jet will always remain in a horizontal 
section. (2) As the section increases in size, the velocity in the 
section becomes less in the inverse ratio. This was a crude form of 
the equation of continuity. From considerations like these Newton 
satisfied himself that the theorem of Torricelli was correct. He 
published his arguments in the Principia in 1687. But the 
theorem was first proved completely by D. Bernoulli in 1738. 
In 1755 Euler gave the equations of motion 

Da 1 dp Dv~ 1 dp d V 

Dt p d . r Dl p dy ' Dl p dz ’ 


where 


D 

Dt 


d 

dt 


dx 



-I- w 


d_ 

dz ’ 


u, v t w being the components of velocity and p the pressure at 
(x, y , z). The equation of continuity was also given by him. 
It was first published in the Memoirs of the Academy of Sciences 
of St. Petersburg in 1770. It is 


dp **(pu) d(pv) d(pw) = 

dt f dx dy dz 
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With these equations Euler proved Torricelli’s theorem in the 
following manner : — 

When the motion is steady, we have - ~ = ^ = 0 . Also, 

J dt St di 

because the only external force is gravity, we have X=Y = 0, 
and Z — g. The equations of motion become 


du du du 1 dp 

U ± V — + W — - = ; 

dx dy dz p dx 


dv dv dv 

U T T- + w *- z 

rlx dy dz 


1 dp 

pdy ; 


dw dw dw 1 dp 

u ^. +v ^ + w fr z =,j -- P r z - 


Since p is constant for ordinary fluids, the equation of con- 

, du dv dw 

tinuity becomes ~ + — -*f — =0. 

dx dy dz 

Further, since u, v, w are the component velocities in an irro- 
tational motion, we have 

d<f> d<b d(j) 

Usz 5T> r = T'i 

dx dy dz 


Hence 


du d l (j> dv 
dy dt/dx dx ' 


And we have two similar equations. 


du dv dw 1 dp . d / _ 1 dp 

*’ - 0r - * 3 ; 


dx 


dx 


Similarly, -l — i/f) ■■ 


• , . d 4 1 d/> 

- — , and \ — ( q 2 ) =g- - — . 

p dp dz p dz 


Multiply these equations by dx, dy, dz respectively and inte- 
grate ; then 

i 9 * = — + gz + C; 

P 

which is a relation between the pressure and velocity. 

Since the hole is small, we can neglect q % at the level of the 
water. This gives 

_ p 

0 = * — ~ + C. 

P 

.*. £ q l = — — 2 4* gz. But — — - = 0, at the hole ; 

P P 

at the hole, \ q 2 -gz, which proves the theorem. 

It will be noticed that the equation of continuity, viz., 

* * 7“ ~ 0 becomes, for irrotational motion, ^ * ~ = 0, 


dx dy dz 


dx 1 dy 2 dz 2 



44 DIFFERENTIAL EQUATIONS IN MATHEMATICAL PHYSICS. 


which is the same as Laplace’s equation about which I spoke 
yesterday. Problems in Hydrodynamics can therefore be treated 
in the same way as problems in Electrostatics or Heat. 

For rotational or vortex motion, the equations of motion 
remain the same as before. We have also 

9A-— . 9 B W <9?; Bu 

By Bz ' ^ Bz Bx * Bx By * 

where 77, £ are the components of the curl. 


Then 

D / £ \ i Bu r) Bu £ Bu 
Dt\p) p Bx + p By f p Bz 

D (v\ £ V £ Bv 

Dt \p) p 8r,+ p By p Bz ’ 

/£\ £ Bw + 7] Bw l Biv 

Dt \ p ) p Bx p By * p Bz 

We can solve these equations in terms of a scalar potential # 
(no longer a velocity potential). 




dx f dif dz\ etc., 


r being 


\/ (z — x') 2 4* (y — y') 2 + (z — z') 1 , 


and 

then 


Bu Bv Bw 
Bx + By + Bz 9 

- ?? 

U Bx h By Bz 9 


— —JJL 

By + Bz Bx ’ 


B0 BG_BF 
Bz * Bx By 


We can thus express u , v, w in terms of F , O , H. The anal- 
ysis is very much the same as that which is used for problems on 
the motion of electricity. This is an important point of view 
and is of special importance in the treatment of viscous fluids. 

The first attempt to study the motion of viscous fluids was 
made by Navier in 1827 and the investigation was extended by 
Poisson in 1831 . But all the investigations which treated the 
problem from the point of view of molecular theory, took the mole- 
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cules as mere points, although they were supposed to be of finite 
size at the start. These investigations were therefore only fruitful 
in creating interest, but lwere not of much scientific importance. 
The first person to dea with the problem with any rigour was 
De St. Venant, who is famous for his work in the theory of 
Elasticity. The equations of motion are 


Du 

Dt 


= X- 



1 dp 
p dx 

1 dp 
p 


v dO 

+ 3 dx* V ' 


V dO 

+ 3 <ty 


+ v V 


where 0 is 



1 dp v dO 
p dz + 3 dz 


+ V V 


the dilatation and equals + v- + ?r • 

dx dy dz 


These equations have a meaning only when the first and the 
secorjd differential coefficients of v, v, ir exist and are finite. The 
continuity of these had also been postulated, in deriving the equa- 
tions. Prof. Oseen showed that the motion can be determined even 

in the case when ~ ~ , — , etc., do not exist at all. He 

dx' dt ' dx 9 

critically examined the equations, just as I have done for the' 
conduction of heat 1 and obtained 

, f f)T( x t r). 

) {r(s,/)-7V,0)}rfr= ^ 

-TT 0 


as the equation of the phenomenon of linear conduction. 

Oseen 2 has obtained certain very interesting results. I read 
out to you a translation of his own views on the subject : “The 
differential equations of Navier admit of being transformed into 
integral equations in which the differential coefficients 

dv dv dw d' 2 u d 2 u dp dp dp 

¥ ’ di ’ ~¥ ’ ¥? ’ ¥}* ’ etc " ¥x'¥ y '¥’ 

do not occur. Indeed it can be shown that the integral equa- 
tions possess systems of solutions for which these differential 
coefficients do not at all exist. One must then ask whether these 


1 *• Constitution of Matter and Analytical Theories of Heat ’* (Abhandbungen 
der Oeselhcha/t der Wissenschaften zu Gottingen, Bd. 2, 190*1). 

* “ Ueber die Bedeutung der Jntegralgleirbungen in der Theorie der Bewe- 
gung einer reibenden, unzusammendriiekbaren Flussigkeit ” (Arlciv for Malematik, 
Astronomi och Fysilc, Bd. 0, 1910). 
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solutions of the integral equations, which are not at the same time 
solutions of the differential equations, possess only a mathematical 
importance or whether some among them have a physical import- 
ance. I undertake the task to show that the latter is the case and 
that therefore not the differential equations of Navier but the 
integral equations should be considered as expressing adequately 
the physical hypotheses which lie at the basis of the theory of the 
motion of a viscous fluid. 

“The direct derivation of the integral equations of Hydro- 
dynamics without the intermediation of Navier’s differential equa- 
tions. appears to me to be of interest from two points of view. 
First, obviously the theory gains thereby in generality, secondly, 
the theory gains much in simplicity. In order to see this, we need 
only consider how complicated the considerations are which will be 
necessary to decide under what circumstances the potential func- 
tions will possess differential coefficients of the second order. 
These questions and others of the same kind occur in Hydro- 
dynamics if one starts from Navier’s differential equations. If one, 
however, obtains directly the integral equations of Hydrodynamics 
all these questions lose their importance. A similar simplifica- 
tion is obtainable in a similar manner in some branches of mathe- 
matical Physics. I may mention as examples the theory of poten- 
tial, and the theory of the conduction of Heat.” 

I give you now the integral equations of Oseen. They bear a 
resemblance to the equations for the vibrating string with an 
infinite number of edges. 

From the principle : 

Amount of motion coming into a closed volume 

= amount of motion due to extraneous forces, 

we get 


j(w, + A , -«,)*>“ j dr^X'(r) (lv + i f/r^ 


X„ (r) ds= 0 


and two similar equations ; 

U v/ tr / du du du\ 

where X —X— ( u — + v - -- + w — l , 

\ Sx dy Oz / 


and X„ = X (B cos nx + X v cos ny + X 3 cos nz. 


Here X w , X y X . are stresses acting in a plane perpendicular 
to the a?-a.xis. In the case of a perfect fluid the stress is parallel 
to the axis. In that case X v and X , are zero. Thus in the case 
of viscous fluid there are nine component stresses to consider. 

The differential equation has thus been replaced by an inte- 
gral equation. There is no great difficulty in the solution of such 
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an equation and we can show that the old solution satisfies this 
equation. 

T now come to Electrodynamics. The equations are 



W 

dz ’ 


4 



d J 

fir ’ 


477 w = 


fi£_fia 
fir By ’ 


da dZ _ fi F 
r/7 fi ?/ fic 


_db_DX^DZ 
dt Dz fir ’ 

__dc_fi] r _fiX 
fi/ fir By 

where ?/, v, ?/> are the components of tlie total current strength at 
(%, y f z), a, /i, y are the components of the magnetic force, a, fi, c 
are the components of magnetic induction and X , Y, Z are the 
components of the electric force. The equations can be solved 
by the use of vectors (F, G , H) and the potential <t> (cf. rotational 
motion in Hydrodynamics) . 

The equations of the theory of electrons are more general, 
and include the equations given above as a particular case. 
Prof. H. A. Lorentz has studied the motion of electrons. The 
equations of motion are 

1 BX 

c, Dt ^ c By dz 5 

1 BY fia fiy 

c Dt + ^ c ~ Dz fir * 

1 DZ w fifi fia 
c Dt ^ c fir By 9 

fi£ fij DZ^ 

fir + By* Dz 

Da DZ BY , f 

etc -> as before ’ 

where p is called the volume density of the electricity and u, v, w 
are regarded as the component velocities of the element of electricity 
which is at the point (r, y, z) at time t. 
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The method of solution is as follows 

<1H dO 

Put 


a 

fly flz 

» 


_r )F dH 


£ 

flz dx 

> 


<1G <1F 


r 

1! 

a>( 
* 1 
1 

3! 


V 

\dF 

dv 

A 

c dt 

dr 

v 

1 <10 

dv 

1 

c fit 

dy 

y 

1 <111 

dV 

/j 

c dt' 

~<lz 


where, F . G, H , are solutions of 


1 <VF u 

r +P " = °> 

C* (9/ c 


and 

1 d l H tv 

V 2 II — + p - = o 

c l dt 1 M c 

also 

-1 

K, 

! 

c^l ^ 
3?| S2 

+ 

II 

© 

and 

df dG dH 1 <1V _ (( 
dx + dy + dz + r dt 


There is no difficulty now in getting a solution. If p~Q f we 
have the equation of wave propagation. From this solution as a 
basis we can build up the complete solution. 

Partial differential equations are found to be quite serviceable 
for all ordinary purposes. But in a rigourous treatment of any 
problem they have to be relegated to a secondary place. The 
work of the last 120 years is of importance as it shows that 
most results and processes dealing with differential coefficients are 
valid only under certain conditions. Thus Taylor’s theorem is not 

1 

always valid. For example, the expansion of e ^” a ) 2 according 
to this theorem, viz.. 
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is wrong because every member of the expansion is zero. This 
discovery, due to Cauchy in 1815, astounded Lagrange who believed 
that if the differential coefficients of all orders existed and were 
finite, the function would be expansible by Taylor’s theorem. But 
T have mentioned to you an example which is in direct contradic- 
tion to his belief. However, in spite of its limitations, the partial 
differential equation will continue to be useful. 

Now I will briefly recapitulate the essence of what 1 have said 
in criticizing the common use of partial differential equations. 
If we have a string with an edge, the first and second differential 
coefficients do not exist at the edge. The equation of motion be- 
comes meaningless in this ease. If we have the initial state of heat 

f(x)~r for #>0, and f(x)~ — x for#<0> 


the first and the second differential coefficients are non-existent at 
the origin. Consequently, the equation of the linear conduction 


of Heat = 
• (it 


<9 2 V 

dj* 


has no meaning for t = 0, x ~ 0. 


In Electrostatics. 


Magnetostatics or (Gravitational Potential, if p = cos log-, v 2 F 

does not exist. In Electrodynamics too a similar remark applies. 
In Hydrodynamics, as Prof. Oseen says, it is not always the 
differential equations but the integral equations which truly express 
the fundamental physical hypotheses. 

It is quite possible that at some time in Hu* future we shall be 
so advanced in oui thoughts that a differential equation will appear 
to us as but a very crude mathematical instrument, and we shall 
discard it in favoui of the more powerful and more refined integral 
equation. 

It has been a real pleasure* to me, gentlemen, to come and 
deliver lectures before you, and 1 shall be much pleased if I find 
that my lectures have been of use in fostering a spirit of research 
here. 









